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Abstract. We present an alternative but intuitive method of obtaining the true or intrinsic quasar luminosity function 
through the observed one in the framework of gravitational lensing. Such a so-called inversion of the magnification 
bias can be straightforwardly performed by means of the Mellin transformation instead of other complex methods. 
Application of our approach to the well-known form of the observed quasar luminosity function identifies the previous 
results in literature. 
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It has been argued for decade whether the magnification bias due to gravitational lensing by the matter clumps in 
the universe affects our quasar number counts N {S), or equivalently the determination of quasar luminosity function 
$ (i, z), and therefore, accounts at least partially for the apparent evolution of quasars (Turner 1980;1981; Canizares 
1981;1982; Avni 1981; Vietri & Ostriker 1983; Schneider 1986;1992; Kayer & Refsdal 1988; Pei 1995). In particular, 
an overdensity of background quasars near foreground galaxies, clusters and even quasars would occur as a result of 
the magnification bias [see Wu (1996) for a recent summary]. All these issues can be simplied as a convolution of the 
true quasar number count N t (S) or luminosity function ${(£, z) around redshift z with a magnification probability 
function p(p) or p(p,\z): 
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where S is the flux threshold of the quasar sample, L, the absolute luminosity and, //, the magnification factor. It is 
evident that the magnification probability distribution should satisfy the following constraints 



/ dp,p(fi\z) = 1, 
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which correspond to, respectively, the normalization and the flux conservation. The question now reduces to how 
to find the true quasar number count N t (S) or luminosity function <& t (L,z) using the observed quantities for a 
given magnification probability distribution p(p) or p(/j,\z), namely, the inversion of eq.(l) or eq.(2). To solve eq.(l), 
Schneider (1992) utilized the Volterra equation of the second kind with a kernel function K(S, x) combined with other 
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mathematical techniques and a specific boundary value, whereas for eq.(2) Pei (1995) performed an expansion of 
the true luminosity function $*(£, z) into the Taylor series coupled with a symbolic operator method to separate the 
variable /i from $t(L, z). These sophisticated methods should be in principle applicable to various matter distributions, 
allowing us to derive the true quasar number count or luminosity function. However, the actual application of these 
methods often turns to be complicated, aside from the unknown details of the magnification probability function. 
Motivated by the importance of eq.(l) or eq.(2) in the study of the associations of angular positions of distant quasars 
with foreground objects (galaxies, groups and clusters of galaxies and quasars), we present an alternative but intuitive 
approach to the inversion of eq.(2). Similarly, this approach can be cquivalently employed for the inversion of eq.(l). 
Our method is based on the Mellin transformation (Titchmarsh 1948). For a give function of f(x), its Mellin 
transformation is an integral of 
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so that the source function itself reads 
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where c should be properly chosen to ensure that the singularities of the function x s /(s) are on the left of the routine. 
The Mellin transformation of eq.(2) is thus 
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in which we have adopted a physically reasonable boundary that the magnification probability function p(p\z) — ► 
for a sufficiently large magnification p. While the right-hand side of the above equation is fortunately a product of 
two Mellin transformations of functions p(p) and <fr t (/z L, z), we have 



$ (s,z) =p(s\z)-$t(s,z), 
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where p(s\z) is the Mellin transformation of p(p\z) and $ t (s, z) and$ (s,z) are the corresponding quantities of $ t (L, z) 
and 3> (L, z), respectively. From eqs.(6) and (8) the true quasar luminosity function can be expressed as 
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As it has been shown, our procedure of inversion of the magnification bias is more straightforward and convenient in 
application than the previous methods. 

To demonstrate how efficiently the present method works, we take the exponential L 1 ' 4 form of the observed quasar 
luminosity function used by Pei (1995) 
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with the luminosity evolution 
L z = L*(l + z)- (1+a *> cxp[-(z - z*f/2al 



(10) 
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here ($*, /3, a, z*, ct») are the parameters fitted by observations. Applying the Mellin transformation of eq.(5) to eq.(10) 
yields 
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where T is the usual T-function with a variable of 4(s — /?). If we adopt the same notation as Pei (1995) by defining 
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then the Mellin transformation of p(fi\z) reads 

p(s\z)=exp[Z(s-l\z)] (14) 

A combination of eq.(12) and eq.(14) gives rise to the Mellin transformation of the true quasar luminosity function 
$ t (s,z) = <t (s, z)\p(s\z)] , and the inverse Mellin transformation of <& t (s,z) results in the true quasar luminosity 
function 

1 /-(/3+l)+ioo $ / L\~ S 

<$> t (L,z) = — 4-4 (— r[4(s-0)]-exp[-Z{s-l\z)]ds, (15) 
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in which we have chosen c = /? + 1. Finally, the ratio of the true quasar luminosity function to the observed one is 
simply 
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-p = cxp(-4m + e m ) x / ds T(4 + 2ms) ■ exp \-2msm - Z(/3 r + 2msf3i\z)] , (16) 
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i.e., the result of Pei (1995) [eq.(39)], where m = ft ln(L/L z ), [3 r = (i, and ft = 1/4. 

Nevertheless, we point out that a quantitative analysis of <3? t /$ depends on the observed quasar luminosity function 
$ (L, z) and the magnification probability function p(fj,\z). Except for some specific forms of & (L, z) [e.g. eq.(10)] and 
p(fi\z), numerical computations should be often employed in finding the Mellin transformations $ (s,z) and p(s\z), 
and hence, the true quasar luminosity function § (L, z). In particular, it is relatively hard to get a simple form of p(/j,\z) 
when the lens exhibits a complicated matter structure (see, for example, Schneider 1992). A detailed investigation 
for various objects as lenses is beyond the scope of this short note. We emphasize that the present method may be 
useful in the study of the associations of background quasars with foreground objects. Recall that the association 
problems, if real, have not been well account for to date in terms of gravitational lensing (Zhu et al. 1996). One of 
the possibilities is to abandon the unaffected background hypothesis, namely, the observed quasar number-magnitude 
relation or luminosity function has probably been contaminated by gravitational lensing according to eq.(l) or eq.(2). 
A further study based on more realistic lensing models, incorporating with the cosmological simulations of formation 
and evolution of large-scale structures, would provide a helpful insight into the problem. 
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